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Abstract. Dynamic properties of molecular motors that fuel their motion by actively 
interacting with underlying molecular tracks are studied theoretically via discrete- 
state stochastic "burnt-bridge" models. The transport of the particles is viewed as 
an effective diffusion along one-dimensional lattices with periodically distributed weak 
links. When an unbiased random walker passes the weak link it can be destroyed 
("burned") with probability p, providing a bias in the motion of the molecular motor. 
A new theoretical approach that allows one to calculate exactly all dynamic properties 
of motor proteins, such as velocity and dispersion, at general conditions is presented. 
It is found that dispersion is a decreasing function of the concentration of bridges, 
while the dependence of dispersion on the burning probability is more complex. Our 
calculations also show a gap in dispersion for very low concentrations of weak links 
which indicates a dynamic phase transition between unbiased and biased diffusion 
regimes. Theoretical findings are supported by Monte Carlo computer simulations. 
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1. Introduction 

Recently much attention lias been devoted to experimental and theoretical studies of 
motor proteins, that are active enzyme molecules that move along linear molecular 
tracks by consuming a chemical energy and transforming it into a mechanical work 
[Ij. Motor proteins play important role in many biological processes [21 |3l H]. While 
conventional molecular motors are powered by the hydrolysis of adenosine triphosphate 
(ATP) or related compounds, it was shown recently that a protein coUagenase utilizes a 
different mechanism for the transport along the collagen fibrils |5l|6]. It fuels its motion 
by a proteolysis, i.e., the collagenase catalyzes the cutting of a protein filament at some 
specific positions, and the enzyme molecule is always found on one side of the proteolysis 
site. Since the collagenase molecule cannot cross the broken site once the cleavage of 
the filament occurred, it leads to the biased diffusion of the motor protein along the 
fibril. 

It was suggested that the unusual dynamics of collagenase can be described by 
the so-called "burnt-bridge model" (BBM) [3 il [H [HI [9]. In BBM the motor protein 
molecule is viewed as an unbiased random walker that hops along the one- dimensional 
lattice composed of strong and weak links. Strong links are not affected when crossed 
by the walker, whereas the weak links ( "bridges" ) might be destroyed ( "burnt" ) with a 
probability < p < 1 when the walker passes over them; and the burnt bridges cannot 
be crossed again. However, theoretical studies of dynamic properties of motor proteins 
in BBM are still very limited. Limiting cases of very low burning probabilities and p — > 1 
case have been analyzed earlier by using spatial continuum approximation [7], which is 
valid only for very low concentrations of bridges. But experimental studies [5l|6] suggest 
that the density of proteolysis sites on collagen is significant. A different approach was 
utilized by Antal and Krapivsky [8j, who investigated discrete-time dynamics of motor 
proteins. They calculated the mean drift velocity of the particle V for periodically 
distributed bridges for the entire range of parameters, while the dispersion D has been 
determined only for the special simple case of p = 1 for both periodic and random 
bridge distributions. However, it was pointed out later [9] that the realistic description 
of motor protein dynamics requires a continuum-time description since motor proteins 
move following Poisson statistics. 

Recently, we developed a theoretical approach to investigate continuous-time 
and discrete-time dynamics of collagen motor proteins in BBM [9] for the periodic 
distribution of weak links. First, it was shown that the special case of j9 = 1 corresponds 
to the motion of the single particle on infinite periodic lattices p2)], for which all dynamic 
properties can be obtained exactly. For the general case of p < 1 we developed a 
simple method by considering a reduced chemical kinetic scheme, and some dynamic 
properties have been explicitly calculated. However, there is a deficiency with this 
method. In this approach the analytical expressions for the dispersion and correlation 
functions have not been obtained, although the simultaneous knowledge of the velocity 
and dispersion is crucial for understanding mechanisms of motor protein transport [HIS]. 
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The goal of the present work is to develop a comprehensive theoretical method that 
allows one to determine all dynamic properties of motor proteins in BBM for general 
sets of parameters. 

2. Model 

In our model the single motor protein is represented by a particle that jumps along an 
infinite one-dimensional lattice with transition rates in both directions assumed to be 
equal to one, as illustrated in Fig. 1. The size of the lattice spacing is also taken to 
be equal to one. The lattice consists of strong and weak bonds. Whereas the random 
walker has no effect on the strong links, it destroys the weak ones with the probability 
p when crossing them. After the bridge is burnt, the particle is always assumed to be 
to the right of it. In addition, it is postulated that all bridges are intact at t = 0. 
Therefore, the trapping of the walker between two burnt bridges cannot occur, the 
particle continuously moves to the right. 



-10 12 N-IN 

Figure 1. Schematic picture for tlie burnt-bridge model of the transport of single 
particles. Thick solid lines describe lattice cells with strong links that cannot be 
burned. Thin solid lines correspond to the weak links that can be burned with the 
probability p when the particle crosses them from left to right. Already burned bridges 
are shown with dotted lines. The particle cam hop with equal rates one step forward 
or backward unless the link is already destroyed. 

The dynamic properties of the molecular motors in BBM strongly depend on the 
details of the bridge burning [9] . There are two possibilities of breaking weak links. In 
the first model, which we term a forward BBM, the bridge burns only when crossed 
from left to right, but it always stays intact if crossed in the opposite direction. In 
the second model, which describes a forward-backward BBM, the bridge is burnt when 
crossed in either direction. Both models are identical for the case of p = 1 because of 
irreversible burning. However, for the general case of p < 1 the walker's dynamics in two 
scenarios, although qualitatively similar, can show some difference [9j. For simplicity, 
in this paper we will only analyze forward BBM, although our method can be easily 
extended for forward-backward BBM. 

Two parameters specify the dynamics of the molecular motors in BBM: the burning 
probability p and the concentration of bridges c. Another important factor that strongly 
influences the dynamic properties is the distribution of bridges [8]. Two possible 
distributions of weak links, periodic and random, have been considered so far [8J. In 
this paper we will investigate the periodic bridge distribution, in which the bridges are 
located at a constant distance N = 1/c lattice spacings apart from each other. Thus 
weak links can be found at the lattice sites with the coordinates kN, where k is integer 
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(see Fig. 1). Periodic distribution of bridges is more realistic for the description of 
dynamics of collagenases O [6] . 

3. Computation of Dynamic Properties 

Our theoretical approach is a generalization of the method proposed by Derrida for 
explicit calculation of dynamic properties of the random walker on infinite periodic one- 
dimensional lattices [10]. It was shown earlier that the original version of Derrida' s 
method can be directly applied to the special case of BBM with p = 1 (with periodic 
bridge distribution) [9] , and it leads to analytical computations of all dynamic properties 
of the system in this special case. However, the direct application of the method cannot 
be done for the case of the bridge burning with the probability p < 1 , since it is unclear 
how the transition rates depend on p. We propose a method that explicitly takes into 
account the presence of the weak links, the feature that is not accounted for in the 
application of the original Derrida's method. 

Let us define PkN+i,y(t) as the probability that at time t the walker is at point 
X = kN + i {i = 0,1, ■ ■ ■ , N — 1) given that the right end of the last burnt bridge is at 
the point yN. We note that y and k > are some integers, and 

+ 00 

RkN+iit)= P{y+k)N+^A^)^ ik>0; t = 0,l,---,N -1), (1) 

y=-OD 

where RkN+i{t) is the probability to find the random walker kN + i sites apart from the 
last burnt bridge at time t. It was found in Ref. [9] that in the stationary-state limit 
(t —>■ oo) the explicit expressions for these functions are given by 



RkN+i{t oo) — RkN+i — RqX^ 



1 - —il-x) 



(2) 



where 



and 



-Ro — / ^ = — ; ^ , (3) 

lp\N - 1)2 + ^pN ^p\\ - c)2 + 4pc 

a; = 1 + ip(iV - 1) - iy'p2(^_i)2 + 4pjv. (4) 
Following the Derrida's method [lOj, we introduce additional auxiliary functions 

+ 00 

SkN+i{t)= iyN+kN+i)PyN+kN+i,y{t), (A; > = 0,l,---,iV-l).(5) 

y=-oo 

The dynamics of the system is governed by a set of Master equations: 
dP (f) 

= PyN^lAt) ' PyNAt) + P E PyN~l,y'{t), (6) 

y'=—oo 



for k = i = 0, and 

dPyN+kN,y(t) 
Jt 



PyN+kN+l,yit) + (1 — p)PyN+kN~l,y{t) — 2PyN+kN,yit) , (7) 
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for k > 1 and i = 0. Also, 

PyN+kN+i+l,y{t) + Py]\f+kN+i-l,y{t) — '2.PyN+kN+i,y{t) (8 



dPyN+kN+i,y{t) 



dt 

for A; > and i = 1, ■ ■ ■ , — 1. It should be noted that substituting the definition ([1]) 
into Eqs. (P) - dS]) leads to the equations for RkN+i that have been already obtained in 
Ref. [9], although through the different route, 

= pY.[RkN-i{t)] + Ri{t) - Ro{t), (9) 
fc=i 

= (1 + R,N+i{t) - 2R,r,{t) (10) 



RkN+i-l(t) + RkN+i+l(t) — '^RkN+iit), (11) 



for A; > 1 and 

dRkN+ijt) 
Jt 

for fc>0 audi = 1,---,A^-1. 

In the similar way, the time evolution of the functions SkN+i{t) can be obtained 
from Eqs. ([5]) and Eqs. ([6]) - ([8]). After some rearrangement it can be shown that 

TO /.\ oo oo 

= S^it) - Soit) +pY: 5a7V-l(t) + P E RaN-lit) - i?i(t), (12) 
a=l a=l 

^^^^ = SkN+i{t)+{'^-p)SkN-i{t)-2SkN{t)+{'^-p)RkN^i{t)-RkN+i{t){13) 

for k > 1 and 

kN+i{ ) _ g^^^.^^(^f'^_^g^^^._^(^f'^_2Si^^_^_^(^t)+RkN+i-i{t)—RkN+i+i{t){l4^^ 



dt 
1,- 

At large times the solutions of Eqs. (fT2|) - (fT^ can be found assuming the following 



for A: > and i = I,- - - , A^- 1. 



form, 

S,{t) = a,t + T„ (15) 

with time- independent coefficients aj and Tj. Substituting Eq. f|T5|) into Eqs. f|T2l) - 
(fTij) yields, 

oo 

ai-ao+pJ2(^aN-i = 0, (16) 

a=l 

OfcAT+i + (1 - p)akN~i - 2akN = (17) 

for k > 1 and 

OfcAf+j+i + o,kN+i-i — 2akN+i = (18) 
for fc>0 audi = 1,---,A^ -1. 



Burnt-Bridge Model 6 

One can easily see that Eqs. (fTBI) - (fTSi) are identical to the system of equations 
for the functions Rj [Eqs. (Q - ffTTl) ] in the stationary-state limit (with dRj/dt = 0). 
Therefore their solutions are the same up to the multiplicative constant, i.e., 

O'kN+i = CRkN+i, (19) 

oo N-1 

where i?fc7v+j is given by Eq. ([2]). Because of the normalization condition J2 RkN+i = 

k=0 i=0 

oo N-1 

1 [9], it follows that C = J2 J2 CLkN+i- To find it explicitly, it is convenient to write the 

k=0 i=0 

equations for Tj which result from the substitution of Eq. ( |T5|) into Eqs. ( |T2l) - dHl), 

oo oo 

ao = Ti-To+pJ2 TaN-l +PJ2 " ^1, (20) 

a=l a=l 

O'kN = Tj^N+l + (1 ~ P)TkN~l — '^TkN + (1 — P)RkN-l — RkN+1, (21) 

for k > 1 and 

(^kN+i = ^fcAT+i+l + TkN+i-1 — '^TkN+i + RkN+i-l — RkN+i+1, (22) 

for > and i = 1, - ■ ■ ,N — 1. Summing over Eqs. (l20l) - (1221) produces 

oo AT-l 

C = X! X! ^kN+i = Ro- (23) 
A:=0 i=Q 



Therefore, according to Eq. f[T9|) 

OfcAf+i = RoRkN+i (24) 

with RkN+i and i?o given by Eqs. ([2]) and ([3]). 

Based on this analysis it is now possible to obtain the expression for the velocity 
of the particle. The mean position of the random walker is the following, 

+00 oo N—l 

{^it))= Y.ll(yN+kN + ^)PyN+kN+^,y{t) = 

?/=— oo k=Q i=0 

oo N~l ( +00 1 oo N-l 

EE E (?/^ + + ^)PyN+kN+^,y{t) = E E SkN+iit). (25) 

fc=0 1=0 [y=-oo J fc=0 1=0 

Therefore the mean velocity is given by 

d d 
V = -rMt)) = i:T.^SkNM{t). (26) 

In the large-time limit it can be shown from Eqs. (|T5|l and (|23|) that 

oo N-l 

V{c,p) = J2J2('kN+^ = Ro, (27) 

k=0 i=0 

with _Ro explicitly given in Eq. ([3]). Note that, as expected, this expression reproduces 
the result for V{c, p) already obtained in Ref. [9] via the reduced chemical kinetic scheme 
method. 
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In order to calculate the diffusion coefficient D{c,p), it is necessary to find functions 
TkN+i from Eqs. fl20l) - fl22l) . To this end, it is convenient to rewrite Eqs. f l20|) - fl22l) by 
replacing OfcAr+j according to Eqs. f l2^ and ([2]), and using the following relations, 

oo r) D 

p 5: i?a7V-l - i?l = ^(1 - X) - i?o, (28) 



(1 — p)RkN-l — RkN+l = RkN+i-1 — RkN+i+1 = ~J^i'^ " x)x'', (29) 

that follow from Eq. ([2]) and Eqs. iQ - ffTTl) in the stationary-state limit. As a result, 
Eqs. fl2Ul) - fl22]) are transformed into the following expressions, 

oo 

Ti-To+pY.T^n-1 + T -Ro = 0, (30) 
TkN+1 + (1 - p)nN-i - 2TkN + Tx'' = (31) 

for k > 1 and 

Tfc^v+^+i + TkN+i-i - 2TkN+i + Tx" + ^(1 - x)Ry = (32) 



for /c > and i = 1, ■ ■ ■ , A^ — 1. In Eqs. (1301) - (13211 we defined a new function F, 

r = ^(l-x)-i?^, (33) 

with Rq and x given by Eqs. ([3]) and (jl]), correspondingly. 

We are looking for the solution of the Eqs. (I30p - (1321) in the following form. 



TkN+i = {Ak+B)x' 



l-^(l-x) 



+ , i(e-l)(l-x 



Tx" — ^ '-R^QX^-Ax^iSA) 



2 6 A^ 

There are three unknown variables in Eq. (I34l) : A, B and A. It is easy to see that 
regardless of their values, Eq. (IMl) automatically solves Eq. ([32D- The i?-term solves 
the homogeneous system [Eqs. (!30!) - (132!) without T- and i?o-terms] as follows from 
the comparison of Eqs. ([30]) - ([MD with Eqs. ([IS]) - ([H]) and Eq. ([2]). As we show 
below, the value of B does not affect the result for the dispersion D{c,p) and it may 
remain undetermined. Therefore, there are two unknowns, A and A, that need to be 
determined from Eqs. ([30]) - ([31]). Substituting Eq. ^ into Eqs. ([30]) and ^ yields 
the system of equations for A and A: 

^{l-p)[l + {N-l)x] + A[-x + {l-p)]= -{l-p)a, (35) 

^^[l + (iV-l)^]77^-AT^ =-^ + Ro, (36) 
N (1 — X)'' 1 — X 1 — X 



with 



c.^^r+ '^-^)f-^) (i^.)< (37) 



while the parameters F, _Ro, x are given by Eqs. ( !33l) . ([3]), and 
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Solving the system of Eqs. (153]) - (15^ results in the following expressions for A and 

A: 

p[l + (A^- l)x](-l + ' ^ 

^ ^ (l-p)[ap + i?o(l-x)^] 

Since i^o, a; (and hence F and a) are functions of (c, p), Eqs. ( l38l) and ( l39l) provide the 
expressions for A(c, p) and A(c, p). Then substituting them into Eq. fl34|) yields the 
formula for TkN+i{c,p). 

It should be pointed out that Eq. (IMI) is not the only possible form of T^N+i that 
solves the system of Eqs. ( |30l) - (1321) . As an alternative, one can use, for example, TkN+i 
which is the same as (1341) except for the last term when Ax^ in Eq. flMl) is replaced by 
iAx^. As a result, corresponding expressions for A{c,p) and A(c, p) can be obtained. 
However, these expressions are much more cumbersome than those given by Eqs. ( l38l) 
and (!39l) . and the final resulting equation for the diffusion coefficient D{c,p), as can be 
shown, does not depend on the specific choice of the functions TkN+i- 

In order to calculate the diffusion coefficient D{c,p) additional auxiliary functions 
are introduced, following Derrida's method [TO] , 

+ 00 

UkN+i{t) = E {yN+kN+iYPyN+kN+i,y{t), k>0, z = 0, 1, ■ ■ ■ , Af-1.(40) 

y=-oo 

The next step is to determine a set of equations that govern the time evolution of Uj(t). 
Corresponding derivation is similar to that of the Eqs. (fT2|) - (fT4l) for functions SkN+i- 
Using the definitions of UkN+i, SkN+i and RkN+i, given by Eqs. (HOi) . ([5]) and together 
with Eqs. ([6]) - ([8]), leads to the following expressions, 

— ^ = f/i(t) - Uoit) - 25i(t) + R^it)+pY^U^N-iit) 

a=l 



+2pY.So,N-i{t)+pY.R^N-iit), (41) 

UkN+l{t) + (1 - p)f/fcJV-l(t) - ^UkNit) - 2SkN+l{t) 

dt 

+2(1 - p)SkN-l{t) + RkN+l{t) + (1 - p)RkN-i{t) (42) 



ci?7fc7v(t) 



for k > 1, and 

'^^A:Af+j(i) 



f^fcAr+i+l(^) + UkN+i-l(t) — '^UkN+iit) — 2SkN+i+l(t) 



dt 

-\-2SkN+i-i(t) + RkN+i+iit) + RkN+i-iit) (43) 

for A: > and i = 1, ■ ■ - , A^ - 1. 

The diffusion constant or dispersion is generally defined as 

d 



D = - lim 

2 t^oo dt 



{x(tY) - {x{t)Y . (44) 
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From the definition fl40p and the summation of Eqs. fl4ip 

d 



it follows that 



oo N-1 



J +00 00 N—1 
y=—oo k=0 i=0 



^ E E u,N-,.{t) = E E 



25o(t) 



i?o(t), 



(45) 



which implies that we do not need to actually solve the system (14T|) - (H3l) to get 
the diffusion coefficient. In the derivation of Eq. (H5l) we also used the fact that 
the probabilities RkN+i{t) are normalized. In the stationary-state limit we have 
So{t) aot + To, where = [from Eq. ^] and Tq = - A [from Eq. ([31])]. Also, 
Ro{t) Rq at steady-state conditions. Therefore, 

(a;(t)2) = 2Rlt + 2B -2A + 2- Rq. (46) 

00 N-l 

Taking into account that in the large-time limit 4(a;(t)) = J2 J2 (^kN+i = Ro [see Eq. 

A:=0 i=0 

f l27|) ] and SkN+i{t) ^ o,kN+it + T^N+i and using the expression for (a;(t)) [Eq. fl25l) ]. we 
obtain 

d d 00 7V-1 

-|(x(t))2 = 2{x{t))-ix{t)) = 2R,{x{t)) = 2i?oE E SkN^t - 00) = 

00 Af-1 00 Af-1 

2-Ro E E WkN+it + TfcAT+i} = 2i?ot + 2i?o E E '^kN+i- (47) 
A:=0 1=0 k=0 i=0 

Substituting the results P6l) and (jUj) into Eq. 
1 



yields the diffusion constant. 



00 N-l 

D = - 2fi - 2A + 2 - i?o - 2i?o E E TkN+i • (48) 

L A,-0 1=0 

Note that in this expression all time-dependent terms disappear, as expected, confirming 
our initial assumptions. 

To obtain the final expression for the dispersion in Eq. (1481) we have to compute 

00 N-l 

the term with J2 J2 T^N+i- It can be shown from Eq. (1511) that 

fc=0 1=0 

00 N-l 00 N-l ( 

EE w. = EE^M(^^ + ^) 



fc=0 i=0 A;=0 i=0 

iil^ ~ I) il - x) 



X 



6 

Utilizing the fact that 

00 N-l 

EE 5^ 

fc=0 1=0 

[as follows from Eqs. ([3]) and 

00 N-l 

E E ^kx 

k=0 i=0 



-R, 







1 (l-x) 



B 



2{l-x 



■[A^ + 1 + {N -l)x\ 



B 

Rq 



and 



l-^(l-x) 



Ax 



2il-x) 



:[iV + l + (iV-l)x] 



(49) 



(50) 



(51) 



Burnt-Bridge Model 



10 



yields Eq. 



taking the form, 

k=0 i=0 



2(1 -x)2 
{N -2){N^ - 1) 



[N+1 + {N- l)x] + 



B 
Ro 



N{N^ - 1) 
6(1 - a;) 



r 



-Rl 



N 



-A, 



24 1-.-' (^2) 

where we also took into account that < x < 1 (for J9 7^ 0) and used the following 

N(N'^-l) ^^'^ 

known results, J2 i{i + I) = — - and J2 — I 

1=0 2=0 

Finally, substituting Eq. (152|) into Eq. 
diffusion coefficient D, 



N{N-2){N^-l) 
4 



2 



-2A + 2 - Ro 



{N -2){N^ - 1) 
24 



Rn 




8]) leads to the following result for the 

N{N^ - 1) 



[A^ + 1 + (iV- l)x] 



6(1 



r 



(53) 



In this expression the undetermined parameter B cancels out and does not affect the 
result for D, as was mentioned above. In the Eq. fl53l) parameters A, A, F, Rq, x are 
given explicitly by Eqs. (l38i) . ( 139|) . ( l33i) . ([3]), and (jl]), correspondingly; the parameter 
a in the expressions for A and A is given by Eq. (!37|) . Since all these parameters are 
functions of only the concentration of bridges c and the burning probability p, Eq. (l53ll 
provides the exact and explicit expression for D{c,p). 



4. Illustrative Examples 

To illustrate our theoretical approach we compute dynamic properties of single motor 
protein in BBM for several sets of parameters. First, let us consider a simple case of 
p = 1, i.e., when any crossing of the bridge burns it. It can be easily shown that 

2 4 2 (N-1) 

^ = °' ^° = ivTT' ^ = NiN + ir "=3(ivTT)- ^''^ 

Parameters A and A have removable singularity aX p = 1 (and therefore so does D). 

Thus we need to compute limA and lim A. To do so we introduce a parameter e -C 1 

p— >i p— >i 

so that p = 1 — e. Then x = j^p^ + 0(6^) as p — > 1. Using this relation we obtain 

, 2iV(2Ar+l) , , 2(iV + 2) 
\imA = — imA = — ) — 55 

P-i 3 (AT + 1)2 ' P-i 3(Ar + l) ^ ^ 

Then, according to Eq. (1531) and Eq. (1271) . the dynamic properties of motor proteins 
are given by 

and 

2 (iV2 + iV + l) _2 (c^ + c+l) 

3 (AT + 1)2 3 (c+l)2 ■ ^ ' 
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These expressions coincide with those obtained in [5] where the velocity and diffusion 
constant were computed for periodic bridge distribution with p = 1 from the original 
Derrida's method, confirming the validity of our theoretical method. 

Another interesting example is c = 1, when every bond in the lattice is a potential 
bridge that can be burned. In this case all auxiliary functions can be easily obtained, 

x=l-^, Ro = Vp^ T = p, a = 0, A=^, A = -^^^tl.(58) 

Substituting these results into Eqs. (l27|l and (153|1 . we obtain for the dynamic properties 
the following expressions, 

V = Vp, D = ^. (59) 

It is interesting to note that in this case the mean velocity depends on the burning 
probability, while the dispersion is independent of p. When p = the diffusion constant 
is equal to one because it corresponds to unbiased motion of the random walker. 
However, any nonzero probability of burning changes the dynamic behavior at large 
times significantly by lowering fluctuations. Thus for any p > the particle is always 
in the biased diffusion regime in the stationary-state limit. 

Generally, in the case of very low burning probability, i.e., when p <^ c, we have 
the following expressions for the auxiliary functions (to the leading order in p/c), 

I (1 - c)p 1 I 1 

yp/c, Rq ^/pc, F ~ pc, a~ — ^ , A rr^^ p j c, A ^ — -.(60) 

Then the dynamic properties are equal to 

V ~ i/pc, D ~ i. (61) 

In the opposite limiting case, when c <C p, our calculations to the leading order in 
cjp produce, 

1 4c^ 2 4 2 

X ~ c( 1), i?o ^ 2c, r ~ — , a ~ -, A ^ -, A ~ — . (62) 

p p 3 3 3 

The corresponding expressions for the mean velocity and dispersion are given by 

\/ ~ 2c, D ~ -. (63) 
3 

The value for the dispersion is rather unexpected since in the case of no bridges (c = 0) 
we have Z) = 1, as expected for the unbiased diffusion of the random walker. 



5. Discussions 



To test our theoretical approach, Monte Carlo computer simulations of dynamic 
properties of motor proteins in the burnt-bridge model have been developed. We 
performed computer simulations by using methods discussed in detail in Ref. [9]. Since 
the velocity of the motor protein in BBM has already been determined in analytical 
calculations and Monte Carlo simulations, we concentrated our efforts on computations 
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of dispersion D{c,p). Diffusion constants or dispersions have been calculated in 
simulations via the formula D = {{x^) - {x)^)/2t. Since the results of computer 
simulations for dispersion indicate large fluctuations, especially for very low c and p, 
typically 10^-10^ simulation runs have been averaged over to decrease stochastic noise. 

The results of Monte Carlo computer simulations for dispersions with different 
parameters are presented in Fig. 2. Comparison of numerical computations with 
analytical predictions indicates a very good agreement for all conditions, and it validates 
our theoretical method. As shown in Fig. 2a, D{c) is a monotonically decreasing 
function for any fixed value of p 7^ 0, and it decreases more slowly for larger values of 
burning probability p. This behavior is expected since the burning of the bridges limits 
the mobility of the hopping particle. At very low concentration of bridges (c <^ p) 
dispersion for all values of p 7^ approached a limiting value of 2/3, in agreement with 
analytical predictions. However, this limiting value is less than D = 1 for the case of the 
diffusion on the lattice without bridges. Note also large fluctuations in D{c) obtained 
from Monte Carlo simulations for c <^ 1 due to the fact that the system, probably, has 
not reached a stationary state yet. 

Dispersion as a function of the burning probability shows a more complex 
dependence, as illustrated in Fig. 2b. For p = the dispersion is equal to one, as 
expected, for any c. The increase in p first significantly lowers the dispersion and then 
D[p) starts to increase. This non-monotonic behavior is rather unexpected because one 
can naively argue that the increase in the burning probability should limit fluctuations 
of the particle, leading to the decrease in the diffusion constant. One can understand 
the dynamics of the single particle in BBM in the following way. Irreversible burning of 
the bridges limits the motion of the particle in one direction, but it allows the particle 
to move further to the right until it crosses another bridge that will burn. It leads to 
larger effective hopping rates, resulting in the larger dispersions. Note that for c = 1 
the dispersions is constant and equal to 1/2 for any p > 0, in agreement with theoretical 
predictions. In addition, D{p) — * 1/2 for very low probabilities of burning {p <^ c) and 

0. 

The presented theoretical analysis provides a full description of dynamics of motor 
protein in BBM. Most observed trends in the dynamic properties can be well understood, 
however there are several surprising observations. First, in the limit of very low 
concentration of bridges (c ^ p) our analytical calculations and computer simulations 
yield D{p) = 2/3 which is not equal to D = 1 expected for the case of c = 0. This 
jump in the dispersion is a signature of the dynamic phase transition [9], that separates 
unbiased diffusion regime (c = 0) from biased diffusion behavior (c > 0). Another 
interesting observation is the fact that D = 1/2 for c = 1 for any value of the burning 
probability p. This can be understood by analyzing the trajectories of the particle. After 
the bridge is burned the random walker spends most of its time near the burning position 
that behaves like a hard wall and it rarely diffuses forward. This lowers significantly 
fluctuations in the system, and the dynamics of the particle become independent of the 
burning probability p. 
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Figure 2. Diffusion constants of the particle in the burnt-bridge model with 
periodically distributed bridges: a) as a function of concentration of bridges; and b) 
as a function of burning probability. Solid lines are results of analytical calculations. 
Symbols are from Monte Carlo computer simulations. 



6. Summary and conclusions 

A comprehensive theoretical approach to describe the dynamics of motor protein 
particles in one-dimensional burnt-bridge model with periodically distributed bridges is 
presented. Our theoretical method allows one to calculate exactly all dynamic properties 
of the system for any set of parameters, and specific calculations are given for the 
mean velocity and dispersion. It is found that dispersion is monotonically decreasing 
function of the concentration of weak links c, while the dependence of D on the burning 
probability p is more complex. Trends in dispersions are discussed using simple physical 
arguments. Theoretical analysis also indicates that there is a gap in the dispersions in 
the limit of c — *• 0, that can be associated with a dynamic phase transition between 
biased and unbiased diffusion regimes. In addition, our calculations show that for any 
concentration of bridges the fluctuations of the particle are lowered. This is due to the 
fact that the particle spends most of the time near the last burned bridge that acts as a 
hard wall. Our theoretical findings are confirmed by extensive Monte Carlo simulations. 

In this paper only periodic distribution of weak links have been considered. This 
situation is realistic for coUagenases moving along collagen fibrils [5] , although different 
distributions of weak links are also generally possible. It will be interesting to investigate 
BBM with random distribution of bridges, for which limited theoretical results are 
available It will also be important to generalize our theoretical approach for studying 
the dynamics of dimers moving on the parallel lattices in BBM. There are theoretical 
predictions that dimers are more effective molecular motors than monomers [6j. In 
addition, our theoretical method provides an explicit analysis of dynamic properties of 
random walkers in quasi-periodic lattice, and this results can be applied for investigation 
of motor protein dynamics in other systems [Ij. 
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